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Abstract 

An upper limit for the r-neutrino mass has been determined from the decay r — ► 5tt v t using 
data collected with the OPAL detector from 1991 to 1995 in e + e~ collisions at yfs ~ M%. A 
limit of 43.2 MeV at 95% CL is obtained using a two-dimensional method in the 57r invariant 
mass and energy distribution from 22 selected events. Combining this result with OPAL's 
previously published measurement using t + t~ — > 31^^ + 3h T i/ r decays, a new combined limit 
of m Vr < 27. Q MeV (95% CL) is obtained. 
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1 Introduction 



The question of whether neutrinos have mass is one of the outstanding issues in particle physics, 
astrophysics, and cosmology. Massive neutrinos are strong candidates for solving the dark 
matter problem of the universe |l|] . Of the three neutrino species, the r-neutrino, is likely to have 
the largest mass. For instance, in the 'see-saw' mechanism |2|] a mass hierarchy exists between 
neutrinos and their corresponding lepton partners, rendering the r-neutrino the heaviest of the 
three known neutrino types. 

On the basis of cosmological arguments a stable r-neutrino with a mass larger than a few 
eV cannot exist 0, however unstable neutrinos may be more massive 0. Previously, OPAL 
has published an upper limit on m VT of 74MeV based on one event in the rare r — > decay 
channel |J in the 1992 data. In this final state the distribution of events in energy E 5n and 
invariant mass of the hadronic system at the two-dimensional limit of the kinematic range 
is sensitive to m UT . 

For this paper, all data collected by OPAL from 1991 to 1995 have been analysed to obtain 
a new limit on the tau-neutrino mass using again the r — > 57r ± z/ T decay channel. Compared to 
the previous analysis, the number of events considered has increased fivefold. 



2 The OPAL detector and simulation 

A detailed description of the OPAL detector can be found in [[J . Subdetectors which are 
particularly relevant to the present analysis are briefly described below. 

The central detector consists of a set of tracking chambers providing charged particle track- 
ing over 96% of the solid angle inside a 0.435 T uniform magnetic field parallel to the beam 
axis. Starting with the innermost components, it consists of a high precision silicon microvertex 
detector, a precision vertex drift chamber, a large volume jet chamber and a set of z-chambers[] 
measuring the track coordinate along the beam direction. 

From 1991 onwards a silicon strip microvertex detector was also present, consisting of two 
concentric layers with readout strips at 50 /im pitch, oriented for azimuthal (0) coordinate 
measurement 0. In 1993 a new silicon strip microvertex detector with z-coordinate readout in 
addition was installed || . 

The jet chamber is designed to combine good space and double track resolution ||, which 
is important for this analysis. It consists of 159 layers of axial anode wires, which are located 
between radii of 255 mm and 1835 mm. The efficiency for separating hits from two adjacent 
particles in the jet chamber is approximately 80% for distances between two hits of 2.5 mm in 
the projection on the r-0 plane @ and drops rapidly for smaller hit distances. The transverse 
momentum resolution of isolated tracks is cr pt /p t = J (0.02) 2 + (0.0015 ■ p t [GeV]) 2 . The jet 
chamber also provides energy loss measurements for particle identification (dE/dx). The dE/dx 
resolution is — 3.2% for minimum ionizing pions in jets with the maximum number 

of hits (159), resulting in a 7r-e separation of at least 2 standard deviations up to momenta of 
14 GeV HQ. 

A lead-glass electromagnetic calorimeter (ECAL) located outside the magnet coil covers the 
full azimuthal range with excellent hermeticity in the polar angle range of | cos#| < 0.82 for 
the barrel region and 0.81 < | cos#| < 0.98 for the endcap region. 

1 The OPAL coordinate system is defined so that z is the coordinate parallel to the beam axis, the radius r 
is the coordinate normal to the beam axis, <f> is the azimuthal angle and 9 is the polar angle with respect to z. 
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The Monte Carlo samples used in this analysis consist of 1.5 million e + e~ — > t + t~ [|T^, 13 
8.5 million e + e~ — > qq" ||14j| and 10 500 e + e~ — > r + T~ff events, which are processed through 
the OPAL detector simulation ||16||. These samples correspond to about 8, 2 and 20 times the 
data luminosity, respectively. 



3 Event selection 

Data collected during the years 1991 to 1995, corresponding to an integrated luminosity of 
155pb _1 and almost 200 000 recorded e + e~ — > r + r~ events have been analysed. The event 
selection is performed in two steps. First, the preselection selects r candidates with five charged 
tracks in a cone. In the second step, background r decays and remaining non-r events are 
rejected. 



3.1 Preselection 



A cone jet algorithm |T7| is employed to assign all tracks and electromagnetic clusters to cones 
with a half opening angle of 35°. For each event exactly two cones are required. The 'signal' 
cone is required to contain exactly five charged tracks with unit total charge. The other 'recoil' 
cone is required to contain at least one track. 

All tracks are required to satisfy the following conditions: 

• Pt > 100 MeV, where p t is the momentum component transverse to the beam direction; 

• at least 20 hits in the central jet chamber. This restricts the acceptance of the detector 
to tracks with | cos^l < 0.963; 

• the distance \d \ of closest approach of the track to the beam axis must be smaller than 
2 cm. The displacement of the track along the beam axis from the nominal interaction 
point at the point of closest approach to the beam must be less than 75 cm; 

• the radial distance from the beam axis of the first hit in the jet chamber associated to a 
track must be smaller than 120 cm. 



To reject non-r events the OPAL standard selection of r pairs is adopted |18|]. The multi 



hadronic background (e + e~ — > qq) is reduced by demanding a maximum of six tracks and 10 
electromagnetic clusters in the event. A cluster is defined as a group of contiguous lead-glass 
blocks which has a minimum energy of 100 MeV in the barrel or 200 MeV in the endcap. The 
requirement on the maximum number of tracks leads to a 5-1 topology of tracks in the signal 
and recoil cone for all preselected events. 



3.2 Final selection 

The background from other r decays and from multihadronic events is reduced by rejecting 
events if the maximum opening angle a max between two tracks in the signal cone is larger than 
10° (see figure |l|a). 

The remaining background is dominated by r decays into three charged particles accompa- 
nied by a photon conversion to an e + e~ pair thus creating a final state with five charged tracks. 
To reject these events the following cuts are applied on the signal side. Events where any track 
has an impact parameter \do\ with respect to the beam axis larger than 0.1cm are rejected 
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(figure 03 ). The minimum transverse momentum p™ in of any track has to be larger than 1 GeV 
(figure 0c). Furthermore the fraction E/p is required to be smaller than 0.7 (figure |l|d), where 
E is the deposited energy in the electromagnetic calorimeter and p is the sum of the momenta 
of the five charged tracks. 

The next two selection cuts exploit the dE/dx information of the jet chamber together with 
the information from the silicon microvertex detector mainly to reject events from r — - > 37r ± 7r°z/ r 
where a photon from the ir° decay has converted. If a track appears to be more likely to originate 
from an electron than from a pion (P e > P n ) or if insufficient dE/dx information is available, 
at least one associated hit in the silicon microvertex detector is required. Here P n (P e ) is the 
X 2 -probability that the track is consistent with the pion (electron) hypothesis derived from 
the dE/dx and momentum measurement (figure|l|e). A test on the total likelihood for the 
57r final state is also performed. The fraction P(57r)/(P(57r) + J2 P(37r)P(e + e~)) must favour 
the 5vr hypothesis (>0.2), where P(5vr) = Ui=i p A i ) and £ P(37r)P(e + e") is the sum of the 
combinatorial possibilities of three particles to be pions and two to be oppositely charged 
electrons (figure 0f). 

Events with a cluster in the electromagnetic calorimeter within the signal cone with energy 
more than 4 GeV and not associated with a track are discarded because this signature probably 
comes from a photon. 

Additional quality cuts on the tracks have been applied to ensure a good reconstruction of 
the 57r system. Each track is required to have at least 40 hits in the central jet chamber. Each 
track fit must have a % 2 per degree of freedom smaller than 2. Given the high density of tracks 
in the 5tc final state, this cut aims to reject events with falsely reconstructed tracks due to 
spatial distortions of the chamber hits or due to hit misassignments by the pattern recognition 
algorithm. Furthermore events are rejected where the angle between a high-momentum track 
(pt > 15 GeV) and any wire plane of the jet chamber is smaller than 0.3°. This cut eliminates 
events with tracks which may be badly reconstructed due to distortions of the drift field in 
direct proximity of the anode and cathode planes. 

After this selection 22 candidate events remain. The positions of these events in the E 5n - 
m 5n plane are shown in figure^. According to the Monte Carlo (MC) simulation, the selection 
efficiency after all cuts is (9.3 ± 0.6)%, where the error is statistical only. 

4 Background 

A high purity data sample is required for an unbiased neutrino mass limit. The background 
can be divided into two classes: (a) r-pair events with a decay misidentified as r — > ^-k^Vt on 
the signal side and (b) non-r events with a topology similar to r decays. The reconstructed 
mass and hadronic energy of these events may be accidentally located close to the kinematic 
boundary in the P 57r -m 57r plane, leading to an artificially low neutrino mass limit. 

The background is estimated from Monte Carlo event samples described in section 0. For 
the t background class (a) we have considered the following decay channels: 

• t — > Stt^Ut-: The r decays into three charged pions one of which undergoes a hadronic 
interaction within the beam pipe or the vertex detector. The final state consists of five 
pions tending to higher invariant masses. 

• t — > Stt^^Ut-: One of the photons from the 7r° decay converts in the detector material 
or a Dalitz decay (tt° — > e + e~7) occurs. If the two electrons are misidentified as pions, 
the reconstructed invariant mass is artificially high. 
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• t — > KgKg7T v T \ If both the Kg decay very close to the interaction point, the final state 
5tt system cannot be distinguished from the signal. The expected bias is small, because 
the mass hypothesis for all tracks is the same as for r — > 57r ± z/ r . 

• t —>■ 57r ± 7r°z/ r : A 20% contamination from these events is expected in the data. As 
explained in section^ these events cannot bias the measurement to lower mass limits and 
this is therefore not a serious background. 

In this background class only one r — > ?>n ± Ti i' T MC event passes the selection corresponding 
to 0.11 data events in the full E^-m^ region. The fraction of m Vr -sensitive or 'effective' r 
background is smaller. An event is denoted as m Ur -sensitive, if its position in the E^-m^ 
plane could lead to a mass limit of below 100 MeV. Based on the E 5n -m 5n distribution of MC 
events in the observed background decay channel it is estimated that less than one tenth of 
these background events would influence the neutrino mass limit. 



Background 
class 


Background 
source 


Expected number 
of events 


Effective number 
of events 


(a) r-> X 


37T ± 7r° 

KgKgV* 37r± 


0.11 ±0.11 
< 0.14(68%CL) 


0.01 ±0.01 


(b) non-r 


qq 
r+r-ff 


0.45 ± 0.45 
< 0.06(68%CL) 


0.04 ±0.04 


total 


0.56 ±0.49 


0.05 ±0.04 



Table 1: Expected background in the selected sample 



Out of the multihadronic MC samples, background class (b), one event is selected. It has 
many clusters in the electromagnetic calorimeter and the reconstructed mass (2.3 GeV) of the 
signal cone is too high to be compatible with a r decay. Normalized to the data luminosity 
this event corresponds to an expected qq background of 0.45 events. For the estimation of the 
effective qq background the multiplicity cuts for tracks on the recoil cone and for clusters in the 
event are relaxed. Then 11 qq MC events are selected. Three of these are located inside the 
kinematically allowed signal region and only one event lies close enough to the boundary such 
that its consideration would have an impact on the extracted limit. It is therefore concluded 
that the qq background that could affect m„ T is only about 0.04 events. 

The expected multihadron background has also been cross-checked using data events looking 
for the 5-2 event topology after relaxing the corresponding multiplicity cut. One such event 
in the data sample is observed with 1.3 expected from the r MC, confirming the direct MC 
prediction of the qq background. 

The background from four fermion events (e + e _ — > r + r~ff ) originates mainly from a r decay 
into three charged tracks combined with the fermion-antifermion pair. In order to estimate the 
contribution of this background, Monte Carlo samples are used where the ff pair is either a qq, 
e + e _ , or /i + /i~pair. No such event passes the selection. 

The total effective background is therefore expected to be about 0.05 events and is considered 
as negligible. The estimated background is summarized in table [I]. 

5 Determination of the mass limit 
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5.1 Likelihood analysis 

The upper limit on the r-neutrino mass is obtained employing a likelihood analysis. The 
probability Pi(rrii, Ei\m v ) for observing each selected event i at the position (mi,Ei) within 
the kinematic plane is derived as function of the neutrino mass m v . To obtain this probability 
the theoretical prediction V(m, E\m„) for measuring the observed distribution in the E^-m^ 
plane is convolved with the experimental resolution R and the detection efficiency e. Hence the 
probability can be written as 

f dm f dEV(m,E\m u )R(m - rrii,E - Ei,a mt ,a Ei ,p)e(m,E) 

Pi(m h Ei\ m v ) - 



J dm J dE V(m, E\ m v )e(m, E) 



The theoretical prediction V(m, E \ m v ) is generated as a function of the neutrino mass m VT 
using KORALZ-TAUOLA (13|, Q including initial-state radiation. The neutrino mass was 
restricted to positive values. The detection efficiency e(m, E) is derived from Monte Carlo 
using full detector simulation |16| . The function used to describe the experimental resolution 
R{m — mi, E — Ei, a mi , a Ei , p) is a two-dimensional Gaussian. The corresponding parameters, 
the errors on the invariant mass a mi , on the energy a Ei , and the correlation p between them, 
are described in the following section. 

The efficiency e(m, E) is, to a good approximation, independent of m and E (e = 0.093 ± 
0.014). Hence the formula for Pi(rrii, Ei\ m„) simplifies to 

Pi(mi,Ei\ m„) = ^ J dm J dEV(m,E\ m u )R(m -m h E- E h er m ., a Ei , p) 

with a constant k — J dm J dEV(m, E\ m v ). 

The sensitivity of the neutrino mass limit to the efficiency is small (see section [7| below). 

5.2 Experimental Resolution 

Most of the events lie well inside the kinematically allowed region in figure ||] such that they do 
not contribute significantly to the mass limit. An event is denoted 'insensitive', if a limit below 
100 MeV cannot be achieved using this event alone. For these events the errors on the track 
parameters are propagated to errors on the invariant mass, energy and the correlation. These 
errors are then taken as input for the binormal resolution function. 

For the other events which are located near or outside the kinematic boundary ('sensitive' 
events), the exact form of the error ellipse in the E^-m^ plane is of crucial importance for the 
determination of the limit on m„ T . Therefore an approach is used which considers the strong 
dependence of the resolution R on the specific topology of the event (i.e., hits in particular 
subdetectors of the tracking system and susceptibility to reconstruction errors). 

For these events the measured four-momenta of each event are used as input for the de- 



tector simulation |TJ| and reprocessed through the full simulation several thousand times. The 
subsample of these events that all have the same reconstruction properties as the original event 
(e.g., same number of tracks with hits in the silicon microvertex detector and in the ^-chambers) 
are used to determine the experimental resolution. This is done by fitting a two-dimensional 
Gaussian function with correlation in an unbinned likelihood fit to the reconstructed E^-m^ 
spectrum of the simulated events. A small non-Gaussian ('tail') fraction of the distributions is 
eliminated to a large extent by discarding events deviating by more than 3 standard deviations 
from the fitted mean. The fraction of discarded events is about 2%, and the remaining events 
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are used to determine the parameters of the resolution function. The typical mass and energy 
resolutions for 5n decays in MC are 20-25 MeV and 500 MeV, respectively. 

In order to possible bias introduced by the tails in the E 57T -m 57T distributions 

the influence of the fraction of events residing in the tails has been estimated. First, a sum 
of two 2-dimensional Gaussian functions has been used in the fit, where the second wider 
Gaussian is introduced to describe the tails. The fraction of this Gaussian has been varied 
by ±50% to estimate the impact of the tails. Alternatively, the sum of a Gaussian and a flat 
pedestal distribution on the E 5w -m 57T plane has been fitted. In both approaches the effect on 
the extracted limit on m Vr does not vary by more than 3.5 MeV. 

It is essential for this analysis that all events taken into consideration for the mass limit 
are well measured. Particularly for decays in which two or more tracks cross within the drift 
chamber volume or approach each other closely, an incorrect hit assignment may cause biases 
in the tails of the E 5n -m 5n distributions. Such biases have been studied using MC events for 
which generated and reconstructed invariant masses and energies can be compared. The same 
procedure as described previously was employed. When the true and reconstructed E 5n -m,5 n 
values do not agree within their errors, striking peculiarities in the distribution after the detector 
simulation are observed. The expected Gaussian peak, generated by the simulation procedure, 
can appear shifted with respect to its input value (figure |3|a) and sometimes ambiguities may 
occur (figure|3|b). 

These effects can largely, but not exclusively, be attributed to pathological track topologies, 
e.g. hit sharing when tracks cross or come very close or when tracks are close to the anode plane 
of the jet chamber. While the shifts as shown in figure|3|a indicate such problematic topologies 
in a clear way, the quality cuts described in section^ are sufficient to remove those pathological 
events in our data sample that are sensitive to vh Vt . From MC simulation it is estimated that 
approximately 50% of the pathological events are rejected by the cut against high-momentum 
tracks in close proximity to a wire plane. 

To assure that this measurement does not deteriorate from such defects, all candidate events 
have been individually inspected. None of the sensitive data events is found to suffer from the 
discussed biases. For the two most sensitive events it has been additionally verified that the 
experimental resolution is nearly constant in the E^-m^ plane. Therefore the resolution has 
been determined for similar events at several positions around the data event (figure |3]c,d). No 
significant deviation from the resolution of the original data event has been found. 



6 Results 

Five sensitive events are retained after the selection described in section [| They are labeled 
by numbers (figure §). Three of them (events 2,3,4) are well reconstructed. All tracks of the 
corresponding signal cones are separated sufficiently and have at least 75% of the maximum 
number of possible hits in the jet chamber. 

Event 1 contains two high-momentum tracks that are close to each other throughout the 
entire volume of the tracking system. Thus the reconstruction of this event is likely to be 
degraded (see position of this event in figure Q). The effect on the limit is small, because the 
event is located in a less sensitive region of the E 5w -m 57T plane. 

In the signal cone of event 5, four track crossings occur in the jet chamber. This results in an 
increased non-Gaussian fraction for the resolution function in the simulated and reconstructed 
events (sect ion |5~2"1) . This fact causes a relatively large error on the invariant mass and energy 



by the likelihood fit. 
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As a result of the background estimation in section^, the probability that one of the re- 
maining five sensitive events is background is 1%. 

The upper limit for m Vr is determined from the 2-dimensional likelihood technique described 



in section |5J] using the combined likelihood of all events which have passed the selection. This 
likelihood function is scanned by changing the assumed true neutrino mass m VT in steps of 
6MeV. A third-order polynomial multiplied by a Gaussian is used to obtain a functional 
description for the likelihood distribution. A 95% CL upper limit of ra„ T < 39.6 MeV is obtained 
by integrating the likelihood function over the physical region of m„ T > 0. The result is shown 
in figure f|a already including systematic uncertainties as described in section [7]. An alternative 
(non-Bayesian) approach using the log-likelihood has been applied as a consistency check, 
resulting in an upper limit in agreement with the one quoted above. 

In table the impact of each sensitive event on the mass limit is listed. 



Event 


Limit variation 
(MeV) 


Limit from this 
event alone (MeV) 


1 


+ 1.8 


97.7 


2 


+7.2 


57.0 


3 


+3.4 


66.6 


4 


+6.1 


58.0 


5 


+0.4 


130.0 



Table 2: Impact of the five sensitive events. The second column shows the effect on the mass 
limit if that event were to be discarded. The last column denotes the limit derived from this 
event alone. 



7 Systematic errors 

The largest systematic uncertainty is from the resolution function, especially from the tail 
fraction. For the sensitive events the parameters of the experimental resolution are varied by 
the average statistical error of the likelihood fit. The effect on the limit is small (0.5 MeV). 



As described in section |5T2| the experimental resolution is determined by the detector sim- 
ulation only for sensitive events. For the remaining events the errors from the track fit are 
used. The consequences of this on and E 5lT are estimated by varying these errors by 30%. 
This is the average deviation observed between the resolution determined by simulation and 
by calculation using track parameter errors. The corresponding variation on the mass limit is 
0.2 MeV. 

The energy calibration has been checked with e + e~ — > events to a level of 5-10 -4 . 

The effect on the mass limit due to this uncertainty is small (0.2 MeV). 

A contamination of 4 events out of the 22 selected events is expected from the decay 
t — ► h'K ±, K v T in the lower E 5n and lower region, as predicted by MC. These events do 
not bias the mass limit to lower values, because the 7r° is not included in the reconstruction. 
To assess a possible impact on the limit, all possible combinations of four events in the lower 
E^- m^ n region have been successively removed from the event sample and the limit has been 
recalculated. The shift in the limit is +0.5 MeV. To be conservative we do not correct for this 
effect. 
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The structure of the r — > 57r ± zy- r decay dynamics only has a very weak effect on the m, 
limit fll9 , [20| , negligible as a contribution to the systematics in this analysis. 



Source 


Limit variation 
(MeV) 


tail fraction description 


3.5 


resolution function for insensitive events 


0.2 


resolution function error 


0.5 


energy calibration 


0.2 


slope of efficiency 


0.3 


tau mass 


0.1 


beam energy 


0.1 


total 


3.6 



Table 3: Systematic effects 



As mentioned in section [| the detection efficiency is assumed to be constant within the 
kinematic region. The effect of a possible m 5n - or .E 57r -dependent efficiency has been taken 
into account by introducing slopes in m 57r and E 5n , varying the efficiency by ±25%. The effect 
on the mass limit was found to be 0.3 MeV. 

A r mass of (1777.0 ± 0.3) MeV has been used [F2T|]; its uncertainty leads to a negligible effect 
on the neutrino mass limit. 

For the beam energy uncertainty, an absolute error of 4 MeV and an energy spread of 28 MeV 



are assumed as obtained by the LEP energy working group [22 



The systematic uncertainties are summarized in tabled All variations of the limit are added 
in quadrature and then added linearly to the mass limit. The likelihood distribution including 
all systematic uncertainties is shown in figure^a; it is obtained by scaling the raw likelihood 
distribution (without systematic effects) in m„ T by the ratio of mass limits obtained with and 
without systematic errors. This likelihood allows us to combine the results of this analysis with 
previous OPAL results on m Vr including systematic errors. 



8 Discussion 

An upper limit for the r-neutrino mass has been derived using the r decay mode r — ► 5-7T v T . 
Including systematic uncertainties, the upper limit 

m Ur < 43.2 MeV is obtained at 95% confidence level. 

This result is based on a data sample that is five times larger than the result previously 
published by OPAL for this channel |5[], and leads to a significant improvement of the limit. 

The combination of this measurement with the previously published OPAL analysis using 
t + t~ — > 3h ± z/ r + 3h T u T decays 0] is obtained by multiplying the respective likelihood curves 
including the systematic uncertainties. The result is shown in figure |]b. From the combined 
new likelihood curve, the upper limit 

m VT < 27.6 MeV is obtained at 95% confidence level. 
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Similar limits of 31 MeV and 30 MeV have been also obtained by the ARGUS p4[ and 



CLEO [^5] experiments, respectively. Recently, a new upper limit of 18 MeV has been deter- 
mined by the ALEPH Collaboration ]BJ| , also by using the results from three- and five-prong 
tau decays. Thus, a tau-neutrino mass of less than 30 MeV is well established and confirmed 
by several experiments. 
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Figure 1: Distributions of the most important quantities used for background suppression in 
the selection. The points with error bars are the data. The histograms denote the Monte Carlo 
expectation, normalized to the luminosity of the data. The cut order and the cut definitions 
are described in section ^T^ . All previous cuts have been applied in each plot. The dashed lines 
indicate the positions of cuts and the arrows point into the selected region. 
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Figure 2: Selected events. The sensitive events are numbered (1-5). The ellipses show the one 
standard deviation contour for the resolution function. The crosses show the errors on and 
E 5n for the insensitive events calculated from the track parameter errors. The kinematically al- 
lowed regions for massless, J^O MeV and 100 MeV r -neutrinos are indicated by the solid, dashed 
and dot-dashed lines, respectively. 
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Figure 3: OPAL simulation of events. Figures (a) and (b) show the output of the simulation 
procedure with one MC event to each plot. The starting point of the arrows shows the true 
position in the E^^m^ plane, the end point the reconstructed position. The latter was used 
as input for the simulation (section 5.2). In cases where a large discrepancy between the true 
and reconstructed values occurs, the simulation shows significant defects such as shifts (a) or 
ambiguities (b) between simulation input and output. Figures (c) and (d) show the analogous 
output for the two most sensitive data events. The crosses denote the input values. The out- 
put distributions are centered around those points and are unambiguous. The stars show the 
additional positions where the resolution was determined. 
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Figure 4: Likelihood versus m Ur for (a) this analysis using r — > ^>-r ± v t and (b) combined with the 
OPAL t + t~ — > 3h ± ^v + 3h T z/ T analysis. The effects of systematic uncertainties are included. 



18 



